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A modification of minimal analytic QCD at low energies 

G. Cvetic a , C. Valenzuela a , and I. Schmidt a * 

a Physics Department, Universidad Tecnica Federico Santa Maria, Valparaiso, Chile 

We present an approach which allows for modifications of the (minimal) analytic QCD model of Shirkov, 
Solovtsov et al. The discontinuity function of the (minimal) analytic QCD coupling parameter is changed at the 
low time-like momenta by a simple ansatz motivated by the vector meson dominance for the vector channel of the 
e + e~ hadronic decay ratio Rv(s). In this way, the additional low energy parameters introduced in the coupling 
parameter make it possible to obtain the correct value of the semihadronic tau decay ratio without having to 
introduce the unusually heavy quark masses of u, d and s quarks, in contrast to the minimal analytic QCD. Some 
numerical results are presented in the leading skeleton approximation for a simple version of such a model. We 
hope that a more systematic analysis of such models leads to an approximately universal coupling parameter. 



00 

o 
in 
o 

a: 

Oh! 



X 



It has been well known that the usual perturba- 
tive QCD coupling a pt (Q 2 ) = a s (Q 2 )/n (where: 



Q 



-<?) 



has nonphysical singu- 



larities outside the time-like momentum region, 
specifically a Landau cut on the space-like axis: 
< Q 2 < A 2 . This can be seen, for example, 
from the expansion 



'pi 



(Q 2 ) 



l 



Ih lnln(Q 2 /A 2 ) 

/3 ln(QVA 2 ) (31 ln 2 (QVA2) 
+G(l/ln 3 (Q 2 /A 2 ) . (1) 



Direct application of the Cauchy theorem gives 
the following dispersion relation for a pt : 



a p t(Q 2 ) = - f 

7T J a 



cr=-A 2 



do pi (a) 
(<J + Q 2 ) 



(2) 



where pi(cr) is the (pQCD) discontinuity function 
of a p t along the cut axis in the complex Q 2 -plane: 
pi(cr) = Ima p t {—a - ie) . 

The physical observables are, by special relativ- 
ity and causality, analytic functions of the asso- 
ciated physical momentum squared q 2 = —Q 2 in 
the entire complex Q 2 plane with the exception of 
the time-like axis (Q 2 < 0). Any QCD observable 
is a function of the invariant coupling a(Q 2 ), so 
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both should have the same analyticity properties. 
Therefore, the cut sector < Q 2 < A 2 in a pt (Q 2 ), 
Eqs. is nonphysical. A simple analytiza- 

tion, which is minimal in a sense, is obtained by 
eliminating the unphysical cut in the dispersion 
relation J5J while keeping elsewhere the same dis- 
continuity function pi(cr), thus leading to a spe- 
cific analytic coupling -Ai(<3 2 ) £0 



*4i(Q 2 ) 



dapiip) 

7T J a=0 {<J + Q 2 ) 



(3) 



Calculational aspects of this prescription were in- 
vestigated in % . Other approaches to analy- 
tization of a p t have focused on the analyticity 
properties of the beta function |8I4| . or the ad- 
dition of power correction terms l/(Q 2 ) n to Ai 
jS] . Our approach considers a modification of the 
pQCD discontinuity function pi (a) of Eq. in 
the a ~ A region of the time-like axis, resulting 
in a different analytic coupling Ai(Q 2 ) 



Ai(Q 2 ) = 



1 



°° da Pl {a) 
=o {* + Q 2 ) 



(4) 



We recall that relation (J2J defines the coupling in 
the entire complex Q 2 -plane except on the time- 
like axis Q 2 < 0. Then, in the analytic approach 
the coupling on the time-like axis s = — Q 2 > 
can be defined in the following convenient way |S] 
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Z7T ' 



(5) 



where the integration is in the Q 2 = cr' plane 
avoiding the (time-like) cut cr' < 0. This rela- 
tion between A\ (Q 2 ) and 21 1 (s) is the same as the 
one between the (vector channel) Adler function 
Dv(Q 2 ) and the e + e~ hadronic scattering cross 
section ratio Rv(s). In addition to relations |i| 
and (0), the following relations hold [7j: 

1 />OC J 

«!(«) = - / -- Pl {v) , (6) 
2 Z" 00 ds&i(«) 



A(Q 2 ) = Q 
"' 2t!(a) = 



o (s + Q 2 )' 



diner 



7T 



(7) 
(8) 



The value of the coupling parameter at the 
origin is in general finite (A(0) = 2li(0) = 

(V 71- ) ir dcr Pi(< 7 )/< 7 -) 

The minimal analytic QCD model 
[cf. Eq. contains only one free parameter - 
the scale A contained in the pQCD discontinuity 
function p\ (cr) . By fitting the measured values of 

high energy QCD observables (at Q or y/s ~ 10 
GeV) to the model, the value A(MS) n/=5 « 
0.26 GeV was obtained 0, corresponding to 
A(MS) n/= 3 ps 0.4 GeV. At lower energies, how- 
ever, a problem appeared. The V-channel of 
the semihadronic r decay ratio, without the 
strangeness production, has been measured with 
high precision by the ALEPH and OPAL jE, Col- 
laborations with the result: i? r .y(A5 = 0) « 1.19, 
and where this quantity has been normalized so 
that in pQCD: i? T: y(AS = 0) = [1 + a pt + 0(a 2 t )]. 
However, fitting this value with the prediction of 
the minimal analytic QCD gave: A(MS) n/ =3 ~ 
0.8 GeV jH]. This problem was then solved by 
the same authors in Ref. 10 . where the value 
A(MS)„.=3 ~ 0.4 GeV was obtained at the price 
of introducing large quark masses m u w wid ~ 
0.25 GeV which implied very strong threshold ef- 
fects in the quantity R T y(AS = 0). The intro- 
duction of (effective) heavy quark masses may be 
regarded as unattractive or hard to justify theo- 
retically. 



A different possibility would be to include 
(some) nonperturbative effects within the time- 
like coupling 2li(s) Ai(Q 2 )]. For that pur- 
pose let us consider first the isovector hadronic 
spectral function i?y(s) [with the normalization 
Rv(s) = l+a p t+0(a 2 )] which has been measured 
with high precision [S]. At low energies \fs ~ 1 
GeV it is dominated by the p-meson resonance 
{M p = 0.776 GeV). On the other hand, the Adler 
function is obtained by the relation 

dsRv{s) 



D V (Q 2 



cr 



(9) 



/o (s + Q 2 ) 2 ' 

which is analogous to Eq. (JJJ . If in the aforemen- 
tioned resonance region we use a narrow width 
approximation (NWA) for Ry(s), in the spirit of 
the Vecton Meson Dominance (VMD), then 

(4^ 2 )- 1 i?y( s )vMD - 2f 2 M 2 S(s-M 2 ) 

+ (47T 2 )- 1 R v (s) pQCI} 9( s -M 2 ut ) , (10) 

where f 2 = (/ 2 )nwa ~ 0.0305, and 6{x) is the 
step function (+1 for x > 0, zero for x < 0). The 
cut value M 2 ut w 1.53 GeV 2 Q3! 

was fixed so that 
the corresponding Adler function J3J is closest 
to the "experimental" Adler function Dy{Q 2 ) exv 
[obtained via relation from the measured 
-Ry(s)ox P ] in the low Q 2 sector (0<Q 2 < 2 GeV 2 ) 
where the quark masses can presumably be ne- 
glected. The agreement of the two Adler func- 
tions is then remarkably close there. Since in the 
analytic QCD R v (s) = l+2Ti(s)+e>((2Ti) 2 ), then 
the VMD form JTHJ) suggests 

M 1 (s) = [R V (s)-l]+0((M 1 ) 2 ) 

» [c f M 2 S(s-M 2 ) - 1] 
+ + 1) 9(s-M 2 ) + 0((2ti) 2 ), (11) 

where 2ti(s) is the minimal analytic QCD time- 
like coupling (containing the parameter A); c/, 
M r and M are additional parameters whose val- 
ues are comparable to 8n 2 f 2 , M p and M cut of 
Eq. I|10|l . respectively. Neglecting higher order 
terms in (|11|) gives 

(s<M 2 ) 
(s>M$) 



2li(s) 



[c f M 2 5(s-M 2 ) 



■1] 



(12) 



and M r < Mo is implicitly assumed. Applying 
the transformation J7J to this VMD-motivated 
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time-like coupling, the corresponding analytic 
space-like coupling is obtained 



-4i(Q 2 ) 



<7 



M 2 Q 2 



Ml 



(Q 2 + M r 2 ) 2 J {Q 2 
1 f°° dapi(a) 



7T } a =Ml (<? + Q 2 ) 



M 2 ) 



(13) 



where 



i 1 f°° d(T t s 
d f = l + - —Pi(cr). 

71" J CT= M n 2 cr 



(14) 



The presented version of the model contains four 
parameters: ct, M r , Mo, and A (the latter is con- 
tained in 2li(s) and in pi(s)). In general, analytic 
coupling l|13(l differs from the minimal analytic 
coupling Ai{Q 2 ), Eq. ©, by terms ~(A 2 /Q 2 ). 

In this presentation we choose, for simplic- 
ity, the additional condition that the two cou- 
plings virtually merge at high momenta, i.e., that 
MQ 2 )~A!(Q 2 ) ~ (A 2 /Q 2 ) 2 . Then, (a) the high 
energy QCD observables (Q, y/s ~ 10 GeV) give 
the value A n .=3 w 0.4 GeV as in the minimal an- 
alytic model [J], and (b) an additional relation is 
obtained 



M 2 c f 



Mr 



+M 2 



dupi (cr) + 



1 [°°da , , , 
Jm% & 



(15) 



Therefore, in such a case the model contains 
two yet undetermined parameters, the scales M r 
and M , i.e., Ai = Ai(Q 2 ;M r ,M ), Mi = 
2ti(s; M r , Mo). These two parameters can then 
be fixed, for example, by the condition that the 
measured semihadronic decay value R T y(AS = 
0) and a "measured" value of the ^-channel Adler 
function D V (Q 2 ) (at a given chosen low Q 2 ) be 
reproduced in the model. In this presentation, we 
will use the approximation of the leading skeleton 
(l.s.), also known as large-/3o resummation. Ac- 



cording to this resummation [12113) 

Dv(Q 2 ) (hs - } = 
1 + 



dtw D (t)a pt {tQ 2 e- 5 / 3 ), (16) 

R r ,v(AS=Q) {u - ) = 

dtW T (t)a^\tm 2 T e- 5 / 3 ),(17) 



1 



where the scheme invariant distribution functions 
twu(t) and fW T (t) are explicitly known |12I13| 
(cf. also ^31) an d are peaked at around t ~ 1. The 
factor e -5 / 3 in the above formulas means that 
the MS scheme is used for et p t, i.e., Aj^g- appears 
there. In Eq. (|17fl . a^^ is a time-like coupling 
|13) . The expressions (|16|l and (|17fl have renor- 
malon ambiguities and the time-like analytic con- 
tinuation ambiguities in pQCD, both manifested 
via the existence of unphysical singularities of a pt 
on the Landau cut. In our model, there is no 
Landau cut, thus no ambiguities, and a p t in the 
above expressions gets replaced by the space-like 
and the time-like couplings, respectively 



D V (Q 2 )^ = 

R T y(AS = 0) { 



dtw D (t)Ai{tQ 2 e^ 5/3 ) 1 (18) 



1 

1+ 4 



^VK T (t)2li(im 2 e- 5/3 ). (19) 



To fix the two parameters M r and Mo, we 
choose as experimental inputs the values 



LV(<3 = 0-7GeV) = 1.00 , 
R T y(AS = 0) = 1.196 



(20) 
(21) 



The value (|20|l was obtained in Ref. ^T] on the 
basis of the ALEPH measurements of Rv(s) [H], 
by applying relation and neglecting the light 
quark masses. The value 1)21(1 is based on the 
measurements of the ALEPH and OPAL Collab- 
orations [5], but with (small) quark mass effects 
subtracted out (cf. also US]). Requiring that the 
leading skeleton expressions 1)181) and 1)19)1 repro- 
duce the values lf2"0"|l and l(2"T]). respectively, gives 
values M r = 0.428 GeV and M a = 0.742 GeV in 
the MS scheme (this corresponds in the V scheme 



4 



G. Cvetic, C. Valenzuela, and I. Schmidt 



to M r = 0.985 GeV and M = 1.708 GeV). With 
these values of the scale parameters, we present 
in Fig. 1 the curve for Dy(Q 2 ) as function of Q 
as predicted by the model, in the leading skele- 
ton approximation (|18fl . and the "experimental" 
curve as obtained by the relation © (in the 
massless quark limit) on the basis of the mea- 
sured values of Ry (s) ■ We see that the two curves 
agree reasonably well in the range Q ~ 1 GeV. 
At higher values of Q, the effects of the c quark 
mass become important |16j . In our calculation 
we used for pi(o~) = Ima pt (— a — ie) the expansion 
JU at the three- loop level [i.e., including terms 
~ln 3 (Q/A).] 



1.2 - 




Q [GeV] 

Fi gure 1 . The "experimental" curve for D V (Q) EH, 
and the prediction of our simple model of VMD-motivated 
analytic QCD in the leading skeleton approximation. 



We see that the model, even in the leading 
skeleton approximation and in one of its simplest 
forms as presented here, reproduces reasonably 
well the low energy behavior of the Adler func- 
tion, and the measured semihadronic r decay ra- 
tio R T , while at the same time the quark masses 
are set equal zero. We recall that the minimal an- 
alytic QCD requires large effective quark masses 
to reproduce the correct value of R T JU] . 

In general, at low energies the quark mass ef- 
fects cannot be entirely neglected. Further, chi- 
ral symmetry breaking, responsible among other 
things for a nonzero pion mass, brings impor- 
tant contributions to some low energy QCD ob- 
servables, such as the axial channel Ra{s) and 



Da(Q 2 )- Such effects will eventually have to be 
taken into account in the framework of the models 
suggested here. Furthermore, the contributions 
beyond the leading skeleton have to be consid- 
ered, possibly by analytization of the power terms 
(apt(Q 2 ))" (n > 2) in the spirit of the approach of 
PP. The type of models suggested here can be im- 
plemented in a more sophisticated way than the 
simplest version presented here - e.g., by going 
beyond the simple VMD-motivated delta ansatz 
for the low energy coupling 2li(s). 
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